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ABSTRACT 

We present photometric and spectroscopic studies of the white dwarf (WD) populations in the 
intermediate-age open clusters NGC 6633 and NGC 7063 as part of the ongoing Lick-Arizona White 
Dwarf Survey (LAWDS). Using wide- field CCD imaging, we locate 41 candidate WDs in the two 
cluster fields: 32 in NGC 6633, and 9 in NGC 7063. Spectroscopic observations confirm 13 of these 
candidates to be bona- fide WDs. We describe in detail our Balmer line fitting technique for deriving 
effective temperatures and surface gravities from optical DA WD spectra and apply the technique to 
the 11 DA WDs in the sample. Of these, only two DA WDs are at the cluster distance moduli, one in 
each cluster. Two more DAs lie 0.75 mag foreground to NGC 6633, raising the possibility that these 
are double degenerate systems in the cluster. If nearly equal-mass binaries, both of these systems 
likely have combined masses above the Chandrasekhar limit. One DB WD is found to be consistent 
with membership in NGC 6633, which would make this the third confirmed He-atmosphere WD in an 
open cluster, though further data are needed to confirm cluster membership. The WD consistent with 
membership in the cluster NGC 7063 has a low mass OAMq), suggesting it may be a He-core WD 
resulting from close binary evolution. Three of the eleven hydrogen-atmosphere WDs in this study 
are observed to have Ca II absorption; the number of DAZs in this study is consistent with previous 
observations that ~ 25% of field WDs are DAZs. 

Subject headings: white dwarfs — open clusters and associations: individual (NGC 6633, NGC 7063) 



1. INTRODUCTION 

White dwarfs (WDs) represent the endpoint of stel- 
lar evolution for the majority of stars. As such, WDs 
present an opportunity to study stellar populations long 
after stellar evolution is complete. WDs in star clusters 
are especially useful due to their residence in a co-eval, 
singie-metallicity population, allowing determination of 
their progenitor masses. 

Studies of WDs in open star clusters are often used 
to address the relation between a star's main-sequence 
mass and that of its WD remnant (the initial-final mass 
relation, or IFMR) . The IFMR represents the integrated 
mass lost by a star over its entire evolution. An un- 
derstanding of this relation is necessary for understand- 
ing chemical enrichment and star formation efficiency in 
galaxies (Fcrrario ct al. 2005) and for understanding the 
origi n and evolution of hot gas in giant elliptical galaxies 
(;e.g.. lMathewill990[ ). The IFMR also plays a vital role in 
the study of the field WD luminosity function (WDLF) , 
as the WDLF represents a convolution of the IFMR with 
the star formation history of a stellar population. The 
WDLF is one of the primary means to determine the age 
of the Galactic thin disk (e.g. JWinget et al.|[T987HWoodl 
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Oswalt et al.lll996l: iLeggett et al.lll998l: Ikilic et all 
yon Hipp el et al. f l2005h . and interest is building in 
attempts to determine the age of the thick disk and t he 
halo via similar methods fe.g.. Ivori Hippel et al. l [200l . 

The first observational work o n WDs in op e n star clus- 
ters was presente d by [R onianisl iiii"fc Angell (|1980D . fol- 
lowed closely by lAnt honv-Twarod (|198ll . I1982 D. Both 
used photographic plates to search for overdensities of 
faint blue objects in intermediate-age open clusters (ages 
~ 100 Myr to ~ 1 Gyr). This work was followed by a 
series of papers by Koester & Reimers, who presented 
follow-up spectroscopy of many of these and other photo- 
graphic WD candidates (Koeste r & Reimers 1981, 1985|, 
I1993L [T996t iReimers fc KoesteJll982l ll988l Il989l [1994) . 
The addition of spectroscopy is crucial. First, it allows 
the unambiguous identification of the bona-fide WDs in 
the samples. In the photometric surveys, candidate WDs 
have been identified spectroscopically as hot horizontal- 
branch, main-sequence or subdwarf stars, and QSOs. 

Second, high quality spectra of the WDs allow the 
physical properties of the WDs to be determined. By 
fitting synthetic spectra to the observed spectra, the 
effective temperature (Tcff) and surface gravity {\ogg) 
of each WD is determined. Using WD evolutionary 
models, these quantities can be converted into a WD 
mass (Mwd) and cooling time (tcooi)- Subtraction of 
Tcooi from a cluster's age results in the total age of 
the progenitor star. From stellar evolutio nary models, 
the pr ogenitor mass is thereby determined. IWeidemannI 
(|2000l) uses data achieved via this method along with 
theoretical stellar evolutionary models to derive a semi- 
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empirical IFMR. 

The advent of wide-field CCD cameras and blue- 
sensitive spectrometers on 8m - f Om class telescopes now 
permits the study of cluster WDs with higher accuracy 
and to large r distances. The first such study is that of 
IClaver et alJ |200I 'l. who use proper motions and digital 
photometry to locate WDs in the core of Praesepe, which 
are then followed up with spectroscopy at large-aperture 
telescopes. 

We have undertaken a similar approach to detecting 
and analyzing WDs in interme diate-age open c lusters, 
with initial results presented in IWilliamj (|2002t l. The 
first results of our study, now tagged with the some- 
what clumsy and outdated (but lAU-approved) appel- 
lation "Lick- Arizona White Dwarf Survey (LAWDS)," 
were se ven massive WDs in t he open cluster NGC 2168 
(M35) (iWilhams et al.ll2004aD . 

At least two other grou ps are conducting co ncurrent 
open cluster WD studies. iKalirai et al.l ()2005bf ) present 
analysis of WDs in the open cluster NGC 2099 (M37) 
and are continuing this work in other older open clusters 
(J. Kalirai 2006, personal communication). In addition, 
iDobbie et all (|2004L [200l l ha ve analyzed s e veral WDs in 
Praesepe not included in the IClaver et al.l (|2001h study. 

In addition to leading to a more defined IFMR, 
the recent cluster observations have revealed the some- 
what surprising possibility that the ratio of hydrogen- 
atmosphere WDs (spectral class DA) to helium- 
atmosphere WDs (the non-DA spectral classes) may dif- 
fer between the field and cl uster WD populations. This 
effect was noted by Reimc rs fc Koester (1994'). and was 
m ore recently explored b y lKalirai et al.l |2005af) . We note 
in IWilliams et al.l (|2006l ) that some of the apparent dis- 
crepancy may be due to the known temperature depen- 
dence of the DA:non-DA ratio, though there is still a 
possible deficiency of non-DA WDs in clusters as com- 
pared to the field. It is therefore important to analyze 
carefully any non-DA spectral type WDs in open cluster 
fields. 

In this paper, we present photometry and spectroscopy 
of WDs in two intermediate-age, sparse open clusters, 
NGC 6633 and NGC 7063. This paper also details our 
spectroscopic analysis technique, wh ich was only briefly 
mentioned in IWilliams et all (|2004a[ ). In §2, we describe 
the two open clusters and our photometric observations 
thereof. §3 describes our spectroscopic observations, 
while §4 presents our spectral fitting routine and test- 
ing of the routine in gory detail. We present our results 
in §5 and discuss these results in respect to the IFMR 
and other open cluster WD issues. In this section we 
also discuss several interesting individual objects. 

In this paper, \ Ye ass ume the solar metallicity value 
of lAsplund et all (|2004 Z = 0.0 13) and the extinction 
curve of iRieke fc Lebofsk^ (|1985f ) with i?y = 3.1. We 
note that initia l analysis of these data were included in 
IWilliamj (j2002D . all analysis presented in this paper su- 
persedes that in the previous work. 

2. PHOTOMETRIC OBSERVATIONS AND 
ANALYSIS 

B — V images of NGC 6633 were obtained from the 
public CFHT archive of the Canadian Astronomy Data 
Centre. These images were t aken as part of the CFHT 
Open Star Cluster Survey (jKalirai et al.l |2001| ) using 



the CFH12K mosaic camera on UT 1999 October 16. 
We reduced the data using the IRAF^ external package 
CFH12K, an extension to the MSCRED mosaic reduc- 
tion package. The twilight flat fields obtained on the 
night of observation were not useful due to high stel- 
lar densities, so twilight flat field images from subse- 
quent nights were retrieved from the archive and used 
to construct flat fields. The data were bias-subtracted, 
trimmed, flat-fielded, and resampled to a common pixel 
sca le and c o ordin ate grid, the basic algorithm suggested 
by IValdesI ()2002f ). Due to the large number of bad 
columns in CCD05 of CFH12K, data on this chip were 
ignored. 

UBV images of NGC 7063 were obtained using the 
Prime Focus Camera (PFCam) on the Lick 3m tele- 
scope on 2001 September 21 and 2002 September 7. PF- 
Cam uses a SITe 2048 x 2048 backside-illuminated CCD 
mounted above a prime focus field corrector and an at- 
mospheric dispersion corrector. The unbinned pixel scale 
is 0'.'296 pixel"^ for a field of view « 10' x 10'. Our ob- 
servations were made with 2x2 binning. 

For the 2001 observations, seeing was ^ 1'.'5, though 
these images suffered from a non-linear response of unde- 
termined origin. The resulting photometry was of suffi- 
cient quality to select objects with UV-excess, but not for 
photometric measurements. For the 2002 observations, 
seeing was « 2'.'2 and the device response was linear, so 
these observations were used for photometric measure- 
ments. 

The log of all observations is given in Tabic [T] 
2.1. Photometric Measurements and Calibration 

Photometry was obtained via p oint-spread fu nction 
(PSF) fitting using DAOPHOT H (lStetsonlll987h . The 
PSF was allowed to vary quadratically in both x- and 
y-coordinates. As the stellar background was crowded, 
stellar PSFs were subtracted from the image, and a sec- 
ond detection iteration was performed. 

Because of shutter timing errors with PFCam and un- 
certainty about the weather during the CFHT imag- 
ing, UBV imaging of portions of both clusters was ob- 
tained with the Nickel 1-m telescope at Lick Observa- 
tory on UT 2002 July 6-8. The Nickel CCD Camera 
has a SITe 1024 x 1024-pixel thinned CCD; pixels were 
binned 2 x 2 for an effective scale of 0'.'56 pix~^ and 
a « 5' X 5' field-of-view. Images were trimmed, bias- 
subtracted, and flat-fielded. We used aperture photome- 
try and the DAOGROW routine to determine total mag- 
nitudes for each star. 

These magnitudes were calibrat ed using a pertu re pho- 
tometry of standard stars from iLandoltj (|1992f) fields 
taken at a variety of times and airmasses. These data 
were used to solve the following transformation equa- 
tions: 

u = U+2.5logU^p+Ao+Ai{U-B)+A2{X-1.25)+A3T (1) 
6 = 5+2.5 log t„xp + Bo + -Bi(-B-V)+B2 (^-1.25) (2) 
« = y+2.51ogt,,p+Co+Ci(B-y)+C2(X-1.25) , (3) 

where u, b, v are the total instrumental magnitudes 
(including a 25 mag zeropoint offset included in the 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 1 
Log of observations. 



UT Date 


Facility 


Instrument 


Seeing 
(") 


Comments 


1999 October 16-17 


CFHT 


CFH12K 


0.7 


Archived images; see Kalirai et al. f2001) 


2001 August 22-23 


Keck I 


LRIS-B 


0.9 - 1.0 


Engineering-grade CCDs 


2001 September 20 


Lick 3m 


PFCam 


1.6 


Non-linear CCD response 


2002 July 6-8 


Lick Im 


Nickel CCD 


2-5 


photometric 


2002 August 6-7 


Keck I 


LRIS-B 


0.6 - 0.7 


variable cirrus 


2002 September 7 


Lick 3m 


PFCam 


2.2 


photometric 


2002 December 8-9 


Keck I 


LRIS-B 


0.9 - 1.2 


clear 


2004 February 12 


Keck I 


LRIS-B 


0.9 


clear 


2005 November 26 


Keck I 


LRIS-B 


0.9 


clear 



DAOphot routines), U, B, V are the standard magni- 
tudes, iexp is the exposure time, X is the airmass, and 
T is the UT time of the observation, in hours. As initial 
T terms for B and V transformations were consistent 
with zero, these are not included in the final calibration. 
Due to a lack of well-measured standards in U, only U 
observations taken on 2002 July 8 could be accurately 
calibrated. The photometric coefficients are given in Ta- 
ble H 

These transformation equations were then applied to 
the star cluster images, and stars with calibrated pho- 
tometry were then used as local standards in the CFHT 
and PFCam imaging to determine the zero points and 
color terms. For the CFHT, we calculated color terms of 
Bi = 0.057 ± 0.029 and Ci = -0.024 ± 0.029. For PF- 
Cam, the corresponding color terms were Ai = 0.095 ± 
0.009, Bi = -0.093 ± 0.005, and Ci = 0.079 ± 0.005. 

2.2. White Dwarf Candidate Selection 

2.3. NGC 6633 

NGC 6633 is a loose open cl uster with recent age es- 
timat es ranging fr om 426 Myr (jPias et al.ll2002l ) to 630 
Myr (jLvngalll987f l. similar in ag e to or slightly younge r 
than the Hyades (« 625 Myr; iPerrvman etaLi I1998D . 
The cluster is slightly metal-poor ([Fe/H] w —0.1), has 
a distance modulus of (to — M)o = 8.01 ± 0.09 (for a 
Pleiades distance modulus of 5.6), and is signi ficantly 
redde ned, with E{B-V) ^ 0.165 ± 0.011 (| Jeffries et aij 
l200^ . 

The color-magnitude diagram for NGC 6633 is shown 
in Fig. [TJ The isochrones shown are interpolated from 
the Z = 0.008 and Z = 0.019 isochrones of lClrardi et all 
(|200l) . assuming scaled solar abundance ratios and mod- 
erate convective overshoot. Interpolation was performed 
using a portion of the StarFISH star-formation history 
code (Harris & Zaritsky 2001, 2004), kindly provided by 
J. Harris. Based on the apparent main-sequence turnoff 
in our new photometry, we estimate a cluster age of 
560tlo Myr {\ogt = 8.75 ± 0.05), in agreement with the 
published values. A clear main sequence can be seen 
down to y « 15 {My « 7). At fainter magnitudes it 
becomes indistinguishable from the very large field pop- 
ulation. 

Figure [1] also shows cooling curves for a range of WD 
masses. As no WD cooling sequence is obvious from 
photometry alone, we use a broad color and photometry 
selection to identify WD candidates. The selection re- 
gion includes DA and DB models at the cluster distance 
and reddening for Mwd > O.4M0 & TcooI < 630 Myr. As 
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Fig. 1. — Color-magnitude diagram for NGC 6633. Grayscale re- 
flects the relative stellar density. Large filled circles are stars in the 
Nickel fields with proper motion membership probabilities > 50% 
in Dias et al. (2001), while small dots are stars in the Nickel fields 
with lower membership probabilities or no measurement. These 
points are also included in the grayscale. Solid isochrones are 
Padova isochrones for logi = 8.7,8.75, and 8.8, interpolated to 
Z = 0.010 and shifted to the cluster distance and reddening. The 
dashed lines indicate candidate WD selection criteria, with the se- 
lected objects indicated as open diamonds. WD cooling curves for 
cluster DAs (solid) and DBs (dashed) are given for logTeooi < 8.8 
and masses of 0.4AfQ(top), O.8M0 (middle), and I.2M0 (bottom). 
No obvious WD cooling sequence is observed. 

we expect any cluster WDs to be younger than the clus- 
ter (and therefore brighter than the faint V limit), these 
photometric selection criteria should include all poten- 
tial cluster WDs. Photometry of all 32 photometric WD 
candidates is given in Table [3l 

2.4. NGC 7063 

NGC 7063 is a poor open c l uster with an age around 
95-125 Myr (i Dias et al.ll2002t iKharchenko et al.ll2005( l. 
making it similar in age to the Pleiades. The intrinsic 
distance modulus i s 9.19 with a f oreground reddening 
E{B-V) = 0.09 (Dias et al.ll2T)0l . We have found no 
metallicity measurements for this cluster, and so make 
the assumption that it is solar. 

The U-B,V B-V,V CMDs for NGC 7063 are 
shown in Figure [21 along with Z = 0.013 isochrones at 
the cluster distance and reddening. The isochrones are 
seen to give a good fit to proper motion cluster members. 
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TABLE 2 

Photometric Transformation Equation Coefficients 

UT Date Filter Zero Point Color Term Airmass Term Obs. Time Term 

2002 Jul 6 B 2.248 ±0.003 -0.083 ± 0.002 0.223 ±0.008 

V 1.986 ±0.005 0.060 ±0.001 0.136 ±0.005 
2002 Jul 7 B 2.611 ±0.016 -0.083 ± 0.002 0.156 ±0.038 

V 2.351 ± 0.009 0.060 ± 0.001 0.070 ± 0.064 

2002 Jul 8 U 4.199 ±0.009 -0.061 ± 0.002 0.511 ±0.012 0.006 ± 0.001 
B 2.573 ± 0.002 -0.083 ± 0.002 0.248 ± 0.008 

V 2.330 ±0.002 0.060 ±0.001 0.145 ±0.008 



TABLE 3 

Candid.we White Dwarfs in NGC 6633 





Object 




RA 




Dec 


V 




uv 


B - 


V 


CTB- 


- V 


Obs. : 


Date 




ID 


Comments 


NGC 


6633 


iLAWDS 


1 


18 


:28: 


:16.7 


6:34:10.8 


17. 


06 


0.02 





,08 





,03 


2002 Aug 


6 


B 


Slitmask 




NGC 


6633 


iLAWDS 


2 


18 


:27: 


:13.6 


6:19:54.5 


17. 


14 


0.02 





,12 





,03 


2002 Aug 


6 


A 


RA 6633-7 


; Slitmask 


NGC 


6633 


iLAWDS 


3 


18 


:27: 


:03.0 


6:24:28.3 


17. 


96 


0.02 





,12 





,03 


2002 Aug 


6 


B 


RA 6633-4 


; Slitmask 


NGC 


6633 


iLAWDS 


4 


18 


:27: 


:10.4 


6:26:15.7 


18. 


82 


0.02 





,16 





,03 


2001 A 


ug 22,23 


DA" 


RA 6633-5 




NGC 


6633 


iLAWDS 


5 


18: 


:26: 


:14.9 


6:29:01.5 


19. 


06 


0.02 





,03 





,03 


2002 Aug 


7 


o 


RA 6633-1 


; Slitmask 


NGC 


6633 


iLAWDS 


6 


18: 


:27: 


:43.7 


6:32:57.1 


19. 


,07 


0.02 


0, 


,38 


0, 


,03 


2002 Aug 


7 


A 


Slitmask 




NGC 


6633: 


iLAWDS 


7 


18: 


:27: 


:49.9 


6:20:51.8 


19. 


,27 


0.02 


0, 


,23 


0, 


,03 


2001 Aug 


23 


DA" 






NGC 


6633: 


iLAWDS 


8 


18: 


:27: 


:23.4 


6:19:49.8 


19, 


74 


0.02 


0, 


,23 


0, 


,03 


2002 Aug 


6 


DAZ 


Slitmask 




NGC 


6633: 


iLAWDS 


9 


18: 


:26: 


:57.5 


6:43:32.0 


20, 


,68 


0.02 


0, 


,15 


0, 


,03 


2001 Aug 


23 


A 






NGC 


6633 


iLAWDS 


10 


18: 


:26: 


:52.4 


6:27:29.3 


20, 


70 


0.02 


0, 


,28 


0, 


,03 


2002 Aug 


6 


B/A 






NGC 


6633: 


iLAWDS 


11 


18: 


:26: 


:29.1 


6:43:22.3 


20, 


,74 


0.02 


0, 


,18 


0, 


,03 


2002 Aug 


6 


B 






NGC 


6633 


iLAWDS 


12 


18: 


:28: 


:49.9 


6:26:12.0 


20, 


79 


0.02 


0, 


,37 


0, 


,03 


2002 Aug 


6 


DA 






NGC 


6633: 


iLAWDS 


13 


18: 


:27: 


:14.9 


6:20:04.1 


21, 


,07 


0.02 


0, 


,14 


0, 


,03 


2001 Aug 


23 


DA? 


Slitmask 




NGC 


6633: 


iLAWDS 


14 


18: 


:27: 


:12.2 


6:21:35.8 


16, 


,93 


0.02 


0, 


,13 


0, 


,03 


2001 Aug 


22 


DB 


RA 6633-6 


; Slitmask 


NGC 


6633: 


:LAWDS 


15 


18: 


:26: 


:08.1 


6:24:51.0 


19, 


,63 


0.02 


0, 


,12 


0, 


,03 


2001 Aug 


23 


DAZ 




































2002 


Aug 


7 




Slitmask 




NGC 


6633: 


iLAWDS 


16 


18: 


:28: 


:04.7 


6:45:06.3 


20, 


,08 


0.02 


0, 


,14 


0, 


,03 


2002 


Aug 


7 


DB 






NGC 


6633: 


iLAWDS 


17 


18: 


:28: 


:03.0 


6:41:27.3 


20, 


,27 


0.12 


0, 


,54 


0, 


,14 














NGC 


6633 


iLAWDS 


18 


18: 


:28: 


:03.1 


6:23:19.4 


20, 


,48 


0.02 


0, 


,51 


0, 


,03 














NGC 


6633 


iLAWDS 


19 


18: 


:27: 


:35.5 


6:38:11.5 


20, 


,55 


0.02 


0, 


,55 


0, 


,03 


2002 


Aug 


7 


QSO 


z = 1.968; 


Slitmask 


NGC 


6633 


iLAWDS 


20 


18: 


:26: 


:56.4 


6:39:05.3 


20, 


,86 


0.23 


0, 


,62 


0, 


,23 














NGC 


6633: 


iLAWDS 


21 


18: 


:27: 


:49.0 


6:45:27.8 


21, 


,02 


0.03 


0, 


,41 


0, 


,04 














NGC 


6633 


iLAWDS 


22 


18: 


:27: 


:07.8 


6:42:29.7 


21, 


,04 


0.02 


0, 


,63 


0, 


,04 














NGC 


6633: 


iLAWDS 


23 


18: 


:26: 


:17.5 


6:22:32.5 


18, 


,77 


0.02 


0, 


,16 


0, 


,03 


2002 


Aug 


7 


O 


RA 6633-2 


; Slitmask 


NGC 


6633 


iLAWDS 


24 


18: 


:27: 


:36.7 


6:26:42.7 


19, 


,02 


0.04 


0, 


,16 


0, 


,06 














NGC 


6633: 


iLAWDS 


25 


18: 


:26: 


:13.6 


6:31:05.5 


19, 


,93 


0.02 


0, 


,27 


0, 


,03 


2003 


Apr 


6 


A 






NGC 


6633 


iLAWDS 


26 


18: 


:27: 


:17.4 


6:22:55.9 


20, 


,09 


0.02 


0, 


,50 


0, 


,03 














NGC 


6633: 


iLAWDS 


27 


18: 


:27: 


:12.3 


6:21:02.1 


20, 


,08 


0.02 


0, 


,29 


0, 


,03 


2002 
2003 


Aug 
Apr 


6 
6 


DAZ 


Slitmask 




NGC 


6633 


iLAWDS 


28 


18: 


:27: 


:22.2 


6:36:26.0 


20, 


,31 


0.02 


0, 


,56 


0, 


,03 














NGC 


6633 


iLAWDS 


29 


18: 


:27: 


:52.6 


6:32:55.0 


20, 


,52 


0.02 


0, 


,44 


0, 


,03 


2002 


Aug 


7 


QSO 


z = 1.643; 


Slitmask 


NGC 


6633 


iLAWDS 


30 


18: 


:27: 


:30.3 


6:24:59.7 


20, 


,78 


0.02 


0, 


,51 


0, 


,03 














NGC 


6633: 


:LAWDS 


31 


18: 


:26: 


:21.6 


6:45:48.0 


20, 


,93 


0.12 


0, 


,39 


0, 


,13 














NGC 


6633: 


:LAWDS 


32 


18: 


:26: 


:16.0 


6:22:58.6 


21, 


,07 


0.03 


0, 


,47 


0, 


,04 


2002 


Aug 


7 


A 


Slitmask 





Note. — Units of right ascension are hours, minutes and seconds, and units of declination are degrees, arcminutes, and arcscconds. 
Coordinates are for Equinox J2000.0. Identifications of B indicate non-WD spectra with He I and H lines; identifications of A indicate 
non-WD spectra with H lines. 

Potential double degenerate star 



and the best cluster age, based mai nly on previously - 
published photometry of HD 203921 (|Hoag et al.lll961l ). 
is 125 Myr (logr = 8.1 ± 0.1). 

Candidate white dwarfs were selected using both 
U — B and B — V criteria, as shown in the color-color dia- 
gram of NGC 7063 (Figure ^ . Comparison of the color- 
color diagram with the color-magnitude diagrams shows 
the utility of multiple colors in selecting WD candidates, 
as many objects that would be selected as WD candidates 
given a single color are eliminated based on two-color 
photometry. In total, nine WD candidates are identified, 
none of which have photometry consistent with a high- 
mass WD at the cluster distance. Photometry of all nine 
candidates is given in Table [H 



3. SPECTROSCOPIC OBSERVATIONS 

Spectroscopic observations were taken between 2001 
August and 2005 November using the blue cha nnel of 
the L RIS spectrograph (Oke et al. 1995; McCarth y et al.l 
119981 ) on Keck I (see Tabic [l]). The 2001 August 
observations used the initial LRIS-B engineering-grade 
2048 X 2048 SITe CCD; all other observations used the 
New Blue Camera, consisting of two 2kx4k Marconi 
CCDs. The 2001 observations therefore have lower sen- 
sitivity in the blue. 

We selected the 400 groves mm" 3400A-blaze grism, 
as it is the available grism with the highest through- 
put for the vital higher-order Balmer lines. The D560 
dichroic was used to permit simultaneous observations of 
the Ha line, though these observations are not presented 
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U-B B-V B-V 



Fig. 2. — Color-magnitude diagrams (left) and color-color plot (right) for NGC 7063. Large filled points indicate stars with proper motion 
membership probabilities > 50% in Bias ct al. (2001). Solid hnes areZ = 0.013 Padova isochr ones of ages l og t = 8.0,8.1, and 8.2. Only the 
logt = 8.1 isochrone is plotted in the color-color diagram. The open circle is HD 203921 from lHoag et aO 1 119611 ). saturated on our images. 
The dashed lines indicate candidate WD selection criteria, with the selected objects indicated as open diamonds. WD cooling curves for 
cluster DAs (solid) and DBs (dashed) are given for logTcooi < 8.2 and masses of O.4M0(top), 0.8Mq (middle), and 1.2Mq (bottom). No 
obvious WD cooling sequence is observed. 



TABLE 4 

CANDm.'iTE White Dwarfs in NGC 7063 



Object 


RA 


Dec 


V 


(TV 


B-V 




U-B 


O'u-B 


Obs. Date 




ID 


Comments 


NGC 7063:LAWDS 1 


21:24:10.3 


+36:26:02.2 


20.03 


0.04 


0.24 


0.06 


-0.34 


0.07 


2002 Aug 


6 


DA 


Slitmask 




















2002 Dec 8 






NGC 7063:LAWDS 2 


21:24:21.5 


+36:26:00.7 


19.20 


0.03 


0.15 


0.05 


-0.48 


0.06 


2002 Aug 


6 


DA 


Slitmask 


NGC 7063:LAWDS 3 


21:24:07.3 


+36:24:45.6 


20.89 


0.06 


0.15 


0.10 


-0.52 


0.11 


2002 Aug 


6 


DA 




NGC 7063:LAWDS 4 


21:24:32.9 


+36:27:51.5 


21.14 


0.04 


0.41 


0.08 


-0.45 


0.09 


2002 Aug 


7 


DA 


Slitmask 


NGC 7063:LAWDS 5 


21:24:43.6 


+36:28:05.6 


22.73 


0.13 


0.43 


0.23 


-0.37 


0.28 


2002 Aug 


7 




Slitmask 


NGC 7063:LAWDS 6 


21:24:49.4 


+36:33:16.4 


21.97 


0.14 


0.42 


0.24 


-0.40 


0.26 


2002 Aug 


6 


DA 




NGC 7063:LAWDS 12 


21:24:33.5 


+36:30:33.5 


22.49 


0.14 


0.39 


0.24 


-0.60 


0.28 


2002 Aug 


7 




Slitmask 


NGC 7063:LAWDS 16 


21:24:50.4 


+36:34:08.2 


21.07 


0.05 


0.48 


0.09 


-0.74 


0.01 


2002 Aug 


6 


QSO 


z = 1.965 


NGC 7063:LAWDS 18 


21:24:20.7 


+36:35:32.0 


20.06 


0.06 


0.27 


0.10 


-0.46 


0.11 


2002 Aug 


6 


A 


Slitmask 




















2002 Dec 8 







Note. — Units of right ascension are hours, minutes and seconds, and units of declination are degrees, arcminutes, and arcseconds. 
Coordinates are for Equinox J2000.0 
^ Insufficient signal to identify spectrum. 



here. Many 2001 observations used a multi-sht mask 
with l"-wide shtlets; these masks were typically not at 
the parallactic angle, and a substantial loss of blue light is 
apparent in these spectra. Subsequent observations used 
a 1" longslit at parallactic angle. The resulting spectral 
resolution (full-width half-max, FWHM) is « 6A. 

We reduced the spectra using the onedspec package in 
IRAF. Overscan regions were used to subtract the am- 
plifier bias, and a normalized flat field was applied to 
the data. Due to very low flux from the flat field lamps, 
the quality of the flat-fielding blueward of w 4000A is 
uncertain. Cosmic rays were removed from the two- 
dimensional spectrum using th e "L. A. Cosmic" Lap lacian 
cosmic ray rejection routine (jvan Dokkuml[200lt ). We 
then co-added multiple exposures of individual objects 
and extracted the one-dimensional spectrum. We ap- 
plied a wavelength solution derived from Hg, Cd and 
Zn lamp spectra. We determined and applied a relative 
flux calibration from longslit spectra of multiple spec- 
trophotometric standard stars. We made no attempt at 



obtaining absolute spectrophotometry for any object. 

Spectroscopic identification of each WD candidate is 
given in Table [3] for NGC 6633 and Table H for NGC 
7063. The major non-WD contaminants in the sample 
are hot subdwarfs (sdB and sdO stars), A- type stars, 
and AGN. We measured the breadth of RS at 20% be- 
low the pseudocontinuum level for a random sub-sample 
of the A-type spectra; these were found to match the 
criterion for field horizontal branch stars (width < 30 A; 
iBeers et"aLll 19881 ). For QSOs, we determined redshifts by 
cross-correlating the spectra with the Sloan Digit al Sky 
Surve y composite QS O spectrum of[y andcn Ber k et al.l 
(|2001h as described in lWilhams et alT(|2004bf ). 

4. MODEL FITTING AND TESTING 

One of the most successful means of measuring DA 
WD properties is the simultaneous fitting of the Balmer 
line profiles in optica l spectra, a method described in 
[Bergeron et al.l fl992D and used in numero u s subsequent 
works (e.g.. [Bergeron et al.l 11994 Il995t iFinlev et all 
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Fig. 3. — Comparison of fit T^ff . Points compare temperatures 
published in Licbcrt ct al. (2005, abscissa) and those from our fits 
to the same spectra, (top) Absolute T^ff values; the dashed line 
indicates equivalence, (bo ttom) Relative Teg differences. Error 
bars are those published in lLiebert et al.l 1)20051 ). The agreement is 
seen to be excellent. 



[19971 iClaverlTalllMltlKalirai et al.ll2005bn . We adopt 
this method for analysis of our WD spectra, though we 
have made some modifications to the aleorithm. The 



major difference is that iBergeron et al.l (jl992f) use the 
Levenberg-Marquardt method to determine the best fit 
and errors, whereas we use a brute-force method, consid- 
ering the entire model atmosphere grid to find the best 
fit and a Monte Carlo simulation to determine the errors. 
We also use a different grid of model atmospheres. These 
and other minor differences necessitate co mparison of the 
solutio ns from our routine and that of IBergeron et al.l 
(|1992f) . In this section, we outline our fitting technique, 
our error analysis, and compare our fits with fits from 
the literature. 

4.1. Spectral Fitting Routine 

The atmospheric models used in the fitting were gra- 
ciously provided by D. Koester, and are slightly mod- 
ified version s of sy nthetic, pure-H atmospheres used in 
iFinlev et al.l ()1997f ). These models cover a range in 
Teff from 10,000 to 80, 000 K and in log 5 from 7 to 
9. These models are interpolated to create a grid with 
ATeff = 100 K and Alogg = 0.02. The models are then 
convolved with a FWHM= gA Gaussian to match our 
spectroscopic instrumental resolution. 

For the line fitting, Balmer-series lines from H/3 
through H9 are considered. For both the model and 
observed spectra, a linear fit is made to the pseudo- 
continuum on either side of the line, and the spectra 




Fig. 4. — Compariso n of fit log g. Points compare gravities pub- 
lished in lLiebert et al.l (2005, abscissa) and those from our fits to 
the same spectra. Very hot (T^ff > 35, 000 K) WDs are shown 
as open symbols; filled symbols indicate WDs with lower tempera- 
tures, (top) Absolute logg values; the dashed line indicates equiv- 
alence, (bot tom) Relative log g differences. Error bars are those 
pubhshed in lLiebert eFalj ||2005|) . For T^g < 35, 000 K, the agree- 
ment is excellent; our fits to hot, low-gravity WDs trend toward 
higher logg values. 

are normalized to the linear fit. The model flux in each 
pixel is determined by averaging the flux at ten equally- 
spaced wavelengths in the pixel; the noise is measured 
empirically from the spectral regions used to fit the con- 
tinuum. The fit of the model to the observations is 
then calculated for each pixel in the one-dimensional ex- 
tracted spectrum. A separate value is calculated for 
each Balmer line i at each T^s and logg in the grid. A 
global minimum Xi min then found for each Balmer 
line. 

For A, (Teff, logg) 



Xf (Teff, logg) 



, the prob- 



ability Pi that a Xi variable exists with a smaller A^ is 



P. = i-r(-,-) 



(4) 



where F is the incomplete gamma functi on and (= 2) is 
the number of degrees of freedom fPress ^ al.lll99^ . In 
other words, Eq.jl] gives the probability that the best-fit 
exists within a contour of constant A^. 

The combined probability 1~P that the best-fit model 
lies outside a given contour of A(Tcff,logg) is then: 



1 -P(Teff, log g) = []l 



(5) 



The global minimum of P{Tcs, logg) then gives the best- 
fitting Teff and logg. 
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TABLE 5 

Comparison of field white dwarf spectral fits 



Object 


TcS o"T„ff logs o"logs T^H^ogg 


Toff, log g 


Tcff.logg 


^cff, logg 


Tcffjogg 




(1) 


(2) 


(3) 


(4) 


(5) 



WU UUUl+4vici 


41800 


90 


8.62 


0.01 










46717, 8.88 




WU 0048+202 


19400 


50 


7.94 


0.01 


20340, 7.97 










20160, 7.99 


WD 0214+5d8 


21000 


50 


7.96 


0.01 


21180, 7.83 








21633, 7.864 




WD Uo4D-Uil 


424UU 


80 


8.98 


0.02 


4U54U, 9.22 






43102, 9.092 








QOfinn 
ozuuu 


f^n 
ou 


O.OU 


U.Ul 


ozyzu, O.OU 


OlDDU 


O. i o 








WD 0352+096 


14800 


50 


8.18 


0.01 




14770 


8.16 








WD 0438+108 


25800 


50 


8.04 


0.01 


27230, 8.05 


27390 


8.07 








WD 0501+524 


61500 


220 


7.58 


0.01 




64100 


7.69 




61193, 7.492 




WD 0549+158 


33300 


70 


7.72 


0.01 








33753, 7.664 


32747, 7.683 




WD 0937+505 


34100 


30 


7.76 


0.01 










35552,7.762 


35830, 7.88 


WD 1105-048 


15700 


50 


7.82 


0.01 








15576, 7.805 






WD 1201-001 


20300 


50 


8.26 


0.01 


19960,8.26 










19770, 8.26 


WD 1451+006 


26000 


50 


7.82 


0.01 


25670, 7.83 






26371,7.926 


26066,7.704 


24930, 7.89 


WD 1737+419 


21900 


50 


7.82 


0.01 










19458, 7.968 




WD 1936+327 


21100 


50 


7.86 


0.01 


21260,7.84 








21240,7.765 




WD 2025+554 


31000 


50 


7.82 


0.01 










30486,7.762 




WD 2046+396 


65700 


270 


7.80 


0.04 










63199,7.766 




WD 2309+105 


58800 


270 


7.80 


0.02 








57990, 8.073 


58701,7.807 


54410, 7.90 


WD 2357+296 


49700 


70 


7.74 


0.01 










49939,7.596 


51960, 7.52 



References. — ( 1) iBergeron et aUlT99a. (2) IBergeron et al.llT995l. (3) IBraeaglia et al.llT995l. (41 IFinlev eFalllTggTl. 
(5') ILiebert et al.ll200a 



In order to reduce the computational time, we use 
an adaptive grid spacing. Starting with a course grid 
spacing, we determine the global best fit and recenter a 
smaller and finer grid on that fit. We iterate this process 
until we reach our model grid resolution of AToff — lOOK 
and A log 5 = 0.02. 

For objects with very high signal-to- noise ratios (S/N> 
240 per resolution element), numerical errors prevent de- 
termination of Toff' and log 5. The reasons for this failure 
are small systematics in flat-fielding and/or continuum 
fitting; due to the very high signal-to-noise, even these 
small deviations result in very large values for which 
P cannot be calculated due to computational numerical 
limitations. This was readily solved by setting a mini- 
mum noise value, such that (7\/ fx^ohs > 0.004. 

To convert from Tcs and log 5 to the WD mass (Mf) 
and cooling age (tcooi), we used evolutionary models 
provided by P. Bergeron. These models include syn- 
thetic photometry and cooling ages for WDs with CO 

cores and thick hydrogen layers [MhIM* = IP 

using WD evolutionary mod els from | Woodl (I1995D for 
Teff > 30, OOOK and from iFontaine etW (|200lD for 
cooler Toff. We only use those models for WDs with 
0.4Mq < MwD < I.2M0; although lower-mass carbon- 
oxygen models are available, any such low-mass WDs are 
likely He-core products of binary star evolution, so use 
of low mass carbon-oxygen models would be improper. 
Mf and TcooI are interpolated from these models using 
the best- fitting Tcff and log 5. We also calculate syn- 
thetic photometric indices for each WD by interpolation 
from the evolutionary models. 

4.2. Spectral Fitting Error Determination 

As the errors in the calculated Toff and log g are corre- 
lated, we choose to determine random errors via Monte 
Carlo methods. By design, such methods also account 
for potential error sources not included in the fitting, 
such as the decreased S/N in the absorption line profiles. 

We begin by convolving the best-fit spectral model 



with the instrumental resolution. The convolved model 
is then multiplied by the spectroscopic instrument re- 
sponse and scaled such that the number of counts in a 
region of the scaled model surrounding H/3 is the square 
of the pixel-to-pixel S /N measured in the same region in 
the observed spectrum. Poisson noise is then added to 
each pixel. Finally, the instrumental response is divided 
out and the noisy model spectrum is fit using our fitting 
routine. 

Nine independent simulations are run for each observed 
white dwarf, and the standard deviations about the mean 
Toff, logg, Mf, and TcooI are calculated; these standard 
deviations are the errors quoted in all fits. The number 
of simulations was selected to maximize the number of 
points used in the statistical analysis while minimizing 
computing time (~ | hr per simulation). iBeers et al.l 
(|1990| ) note that, for sample sizes similar to our nine 
simulations, the standard deviation does not always give 
a good estimate of the sc ale of the scatter; they find that 
the "gapper" method of IWainer fc ThissenI (l976f ) best 
recovers the true scatter. We calculate the scatter of 
Toff and log g using the gapper method and find identical 
values to within our quoted precision. We therefore claim 
that the nine simulations are sufhcient in recovering the 
measurement errors of Toff and logg. As the smallest 
model spectral grid has a step size of ATeff = 100 K and 
Alogg = 0.02, we adopt minimum errors of (XTeff = 50 
K and criogg = 0.01. 

4.3. Comparison to Published WD Atmospheric 
Parameters 

In order to further test our Balmer line fitting proce- 
dures and derive atmospheric parameters for WDs, we 
make two comparisons of our results with those in the 
literature. First, we fit spectra of WDs from th e Palo - 
mar Green (PC) Survey study of iLiebert erall (|2005D 
provided by J. Liebert. This provides a direct compar- 
ison of the output fr om our fitting routine and that of 
IBergeron et al.l ()1992f ) . This fitting also provides a direct 
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comparison of the at mospheric models of P . Bergeron and 
collaborators used in lLiebert et alJ ()2005l ) and the atmo- 
spheric models of D. Koester and collaborators used in 
our fitting. 

Rather than attempt to fit all 348 DA WDs in the PG 
sample, we selected a sub-sample of 50 objects, chosen to 
cover the range of our model grid in Tes and logg. This 
sample was then pared to exclude three known magnetic 
WDs and two WDs with composite spectra. Two hot 
WDs were also removed due to the presence of He II 
lines, indicating these are DAOs. After these cuts, 44 
objects remained in the sample. 

Figure Bl compare s the T^g from our fits with those of 
iLiebert et al.l (|2005D . With few exceptions, agreement is 
excellent. At high temperatures (Toff > 60,000K), our 
fit temperatures appear systematically higher than the 
published values. This systematic is not of great concern, 
as these temperatures are known to be biased due to the 
prese nce of significant metal opacity in many hot WDs 
(e.g.. iWolff et al.lll998[ ). and because these objects cool 
so rapidly that even a large fractional error in TcooI leads 
to a small fractional error in the progenitor star's age, 
resulting in only small errors in the derived initial mass. 
The most discrepant T^ff measurement at cooler T^s is 
for PGI255-f426 {AT^s = 6800K; ctt„„ = 1135K); our fit 
is poor due to the low S /N of the observation. Excluding 
the hot WDs and PG1255-I-426, the standard deviation 
in the offsets between the published T^s and our fit 
Teff isa(ATeff) = 685K. 

A comparison of the best-fitting surface gravities (Fig- 
ure 2]) reveals a larger scatter than is seen in the temper- 
atures [a{Alogg) = 0.17], though the overall agreement 
is good. There is some evidence for a systematic discrep- 
ancy for hot WDs with log 5 < 7.75. This is due to the 
weakness at high Tcff of the higher-order Balmer lines, 
which provide the most leverage for surface gravity de- 
terminations. Excluding WDs with Tcff > 35000 K, the 
standard deviation in the difference of the log g determi- 
nations is 0.08. 

In summary, our fitting method and model atmo- 
spheres recover very similar atmospheri c parameters 
as the method and atmospheres used in ILiebert et al.l 
()2005l ). with the exception of the lower gravity (log 5 < 
7.75), hot (Tefi > 35,000K) WDs. This agreement gives 
us confidence that our fitting routine is working as ex- 
pected. 

We now compare fits of bright WD spectra obtained 
during our spectroscopic runs with published atmo- 
spheric parameters of these well-studied objects. This 
comparison identifies potential systematic effects intro- 
duced by differing instruments and data reduction tech- 
niques. 

We obtained spectra of 19 WDs, each with atmo- 
sph eric parameters published in one or more of five stud - 
ies: iBergeron et all (ll992lll995D:lBragaglia et all (I1995D: 
iFinlev et all (|1997D : lLie'bert et all (|2005D . Our fits and 
those from the literature are given in Table O We 
noted no obvious systematic in temperature determina- 
tions as compared with previous studies, though the scat- 
ter [(T(ATcff) = HOOK] is significantly larger than the 
quoted internal errors. There is a systematic offset in 
surface gravity for logg > 8.5, with our fits having sig- 
nificantly lower values than previous studies. However, 
only a few high-gravity WDs were investigated, and two 
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Fig. 5.— Spectra of WDs in the field of NGC 6633 and NGC 
7063. The top two spectra are DB WDs; the bottom WD spec- 
trum was unable to be fit due to poor sky subtraction in the Balmer 
line cores. All other spectra are arranged in order of increasing fit 
Tcff (bottom to top). Spectra are normalized to unity at 5300A and 
offset by arbitrary amounts. No correction for interstellar red- 
dening has been applied. Labels are in the format NGC:ID; e.g., 
6633:16 = NGC 6633:LAWDS 16. 

of these WDs have published logg > 9.0, the highest 
gravity in our model grid. The scatter in the logg de- 
terminations is found to be a{A\ogg) = 0.12. 

5. DISCUSSION 
5.1. Candidate White Dwarf Spectral Fits 

Spectra of the WDs in both cluster fields are shown in 
Figure ini and each object is also identified in the cluster 
CMD (Figure Using our spectral fitting technique, 
we fit the spectra of all DA WDs detected with LRIS. 
We present the results of these fits in Table |6| and show 
the Balmer line profiles in Figure [71 We then inte r polate 
using the evolutionary models of Fontaine et al. (2001) 
and photometric models of Holbcrg & Bergeron (200(|) 
to determine WD masses and cooling ages. We note a 
persistent error in continuum fitting of the H9 line in 
NGC 6633:LAWDS 8 (visible in Figure HI); the H9 hue is 
excluded from that fit. 

Cluster membership is determined based on the ap- 
parent distance modulus (m — M)v and cooling age of 
each WD. Potential cluster members are defined as those 
with Tcooi less than la above the adopted cluster age 
and distance moduli within 2a of the cluster distance. 
For NGC 6633, these criteria are TcooI < 625 Myr and 
8.34 < (m — M)v < 8.70. In the literature, no errors 
are explicitly stated on the distance to NGC 7063, so we 
adopt the scatter in published values of ^ 0.2 mag; the 
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TABLE 6 
DA White Dwarf Spectral Fits 



Object 


Dcs 


gnation^ 


S/N^- 


(K) 


(K) 


log 9 


"logs 


Mq 


Mq 


^eool 
log(yr) 


--eool ( 

log(yr) 


m - M)y 




NGC 6633:LAWDS 4 


WD 


,11827+064 


114 


20800 


50 


8.26 


0.01 


0.79 


0.01 


8.035 


0.010 


7.75 


0.02 


NGC 6633:LAWDS 7 


WD 


J1827+063.1 


114 


18500 


50 


8.40 


0.02 


0.87 


0.01 


8.314 


0.013 


7.78 


0.03 


NGC 6633:LAWDS 8 


WD 


,11827+063.2 


62 


11400 


50 


8.84 


0.03 


1.11 


0.02 


9.255 


0.012 


6.55 


0.06 


NGC 6633:LAWDS 12 


WD 


,11828 + 063 


97 


17500 


90 


8.12 


0.02 


0.69 


0.01 


8.168 


0.017 


9.64 


0.04 


NGC 6633:LAWDS 13 


WD 


J1827+063.3 


34 


50000 


1110 


7.24 


0.16 


0.43 


0.04 


5.674 


0.388 


13.29 


0.31 


NGC 6633;LAWDS 15 


WD 


,11826+064 


167 


27200 


50 


7.72 


0.01 


0.49 


0.01 


7.098 


0.003 


9.90 


0.02 


NGC 6633:LAWDS 27 


WD 


J1827+063.4 


79 


16100 


90 


8.24 


0.02 


0.77 


0.01 


8.370 


0.015 


8.61 


0.04 


NGC 7063:LAWDS 1 


WD 


J2124 + 363.1 


113 


15500 


60 


7.50 


0.01 


0.37 


0.01 


7.930 


0.008 


9.61 


0.04 


NGC 7063:LAWDS 2 


WD 


J2124 + 363.2 


162 


15300 


50 


7.86 


0.01 


0.53 


0.01 


8.169 


0.008 


8.20 


0.03 


NGC 7063;LAWDS 3 


WD 


J2124 + 364.1 


111 


17900 


140 


8.04 


0.03 


0.65 


0.02 


8.076 


0.028 


9.90 


0.08 


NGC 7063:LAWDS 4 


WD 


J2124 + 364.2 


22 


18300 


190 


8.80 


0.06 


1.09 


0.03 


8.629 


0.053 


8.90 


0.13 
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Fig. 6. — Color-magnitude diagram for NGC 6633 with spec- 
troscopic identifications. Grayscale reflects the relative stellar den- 
sity in the CFHT imaging, while points are all stars in the Nickel 
fields. The logt = 8.75 isochrone is shown (solid), along with cool- 
ing curves for for cluster DAs (solid) and DBs (dashed) are given 
for logTcooi < 8.8 and masses of O.4M0(top), 0.8Mq (middle), and 
1.2M0 (bottom). The photometric selection criteria are dashed 
lines. Crosses indicate objects identified as non-WD or objects 
not spectroscopically identified. Small circles indicate non-cluster 
member DAs (filled) and DBs (open). The large filled circle is the 
cluster DA, while the large open circle is the potential cluster DB. 
The open squares are the two potential double degenerates. 



corresponding selection criteria for NGC 7063 WDs are 
Tcooi < 158 Myr and 9.07 < (m - M)v < 9.87. Apply- 
ing these criteria, we find one cluster member in NGC 
6633 (NGC 6633:LAWD S 27), thou gh up t o three more 
members may exist; see i i5.2.1l and tj5.2.2l One cluster 
member in NGC 7063 (NGC 7063:LAWDS 1) is detected. 
A graphical exposition of this process is shown in Figures 
El and El 

We calculate progenitor masses for cluster mem- 
ber WDs using stellar evol utionary track s (incl uding 
convective overshoo t) from iGirardi et al.1 (|2000f ) and 
iBertelli et al.l (|1994f ). The calculated TcooI is subtracted 
from the cluster age to produce the progenitor star's life- 
time from the zero-age main sequence (ZAMS) through 
the planetary nebula phase. The ZAMS mass corre- 
sponding to the progenitor lifetime is then determined 
from the evolutionary tracks. As the cluster metallicities 
differ from those of the evolutionary tracks, we linearly 
interpolate the progenitor masses to the assumed metal- 
licity. Progenitor masses for the cluster WDs are given 



in Table [71 

There may be some field WDs that meet our selec- 
tion criteria for cluster members. We estimate the ex- 
pected number of interlopers as follows. The bright WD 
luminosity fu n ction from the Palomar-Green sample of 
iLiebert et all (|2005l) can be approximated by a power 
law: 

log = 0.84A/y - 13.3, (6) 

where (f> is in units of pc^^ 0.5 mag~^. We integrate 
this to limiting magnitudes of My = 12.5 for NGC 6633 
and Mv = 11.75 for NGC 7063, both of these values 
approximating the luminosity of a 1.2 A/© WD with a 
cooling age equal to the cluster age. The resulting field 
WD densities are 1.6 x lO'^ pc'^ for NGC 6633 and 3.9 x 
10"^ pc-3 for NGC 7063. The volumes in which these 
field white dwarfs could reside are truncated pyramids 
with faces defined by the image boundaries and bases 
defined by the distance selection criteria for each cluster. 
Assuming all extinction is foreground to the selection 
volume, these volumes are 1060 pc^ for NGC 6633 and 
520 pc'^ for NGC 7063. So, we expect an average of 1.7 
field WDs in the NGC 6633 and 0.2 field WDs in the 
NGC 7063 to meet the cluster member selection criteria. 

Assuming Poisson statistics, there is an « 18% proba- 
bility that the NGC 6633 region has no field WDs, while 
there is an « 82% probabifity that NGC 7063 has no field 
WDs. Based on these relatively ambiguous statistics, we 
cannot claim with any confidence whether the candidate 
cluster WDs are truly cluster members or are just part of 
the field WD population. Proper motion measurements 
will be necessary to confirm either scenario. 

5.2. Notes on Individual Objects 
5.2.1. NGC 6633:LAWDS 4 and NGC 6633:LAWDS 7 

These WDs are both relatively massive (O.79A/0 and 
0.87AfQ, respectively), as one might expect for NGC 
6633 WDs such as NGC 6633:LAWDS 27 (0.771/0). 
These two WDs also have nearly identical distances 
[(to - M)v = 7.75 ± 0.02 and (m - M)y = 7.78 ± 0.03], 
placing them foreground to the cluster by 0.77 and 0.74 
mag, respectively. These two WDs therefore make strong 
candidates for binary WDs. 

From our data, it is not yet possible to tell if one 
or b oth of these ob j ects i s actually a double degener- 
ate. iBergeron et al.l ()1989l ) find that the unresolved bi- 
naries can be fit acceptably by a single-star spectrum, 
with the resulting fit logg and Toff intermediate to that 
of the two binary components. The color of NGC 
6633:LAWDS 4 {B-V = 0.16 ± 0.03) is what would 
be expected based on the spectral parameters and the 
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Fig. 7. — Balmer line profiles for all DA WDs in the cluster fields. The plotted lines range from H/3 (bottom) to H9 (top), normalized to 
the pseudo-continuum just outside the plotted regions and arbitrarily offset vertically. Thick solid lines are the best-fitting models, while 
the thinner histogram is the observed spectrum. Labels are in the format NGC:ID; e.g. 6633:4 = NGC 6633:LAWDS 4. 



(Reddened) photometric models of iHolberg fc BergeronI 
(I2OO6,) : {B-V)modci = 0.11, but the color of NGC 6633: 
LAWDS 6{B-V ^ 0.38 ±0.03) is redder than the mod- 
els would imply — model = 0.17]. This redder color 
could suggest a cooler component to the system. 

The presence of binary WDs in open clusters has 
been predicted bas ed on dynamical models (e.g., 
iHurlev fc Sharall20f)3h . Further, the presence of such bi- 
naries may be useful for explaining the perceived deficit 
of WDs in some open clusters and for determini ng the 
distribution of binary mass ratios (Williams 2004) . 

However, if these WDs are double-degenerate cluster 
members, the binary components would have to have 
virtually identical luminosities to appear 0.75 mag over- 
luminous, further implying that the two WDs are likely 
of similar masses. Since WDs of different masses cool at 
different rates, it would be a surprising coincidence if two 
separate systems both happened to have different-mass 
components presently with the same luminosity. 

Further, if these are nearly equal-mass component bi- 
naries, both pairs would have a combined mass above the 
Chandrasekhar mass, making these candidates for Type 
la supernova progenitors, depending on their orbital sep- 
arations. However, the ESO Supernova la Progenitor 
SurveY (SPY) has only detected 3 binaries with com- 
bined masses over the Chandrasekh ar limit out of over 
100 detected binary WDs (Napiwo tzki et al.|[2005.) . We 
therefore emphasize that the identification of these ob- 



jects as double degenerates is, at present, speculative. 

It may be possible to determine whether these two ob- 
jects are binaries with additional, high-resolution spec- 
troscopy of the Ha line, where two non-LTE line cores or 
a single core with variable velocity should be detectable. 
In addition, second-epoch deep imaging of this cluster 
will allow proper motions to be measured and the prob- 
ability of cluster membership to be determined. 

5.2.2. NGC 6633. -LAWDS 16 

The spectrum of NGC 6633:LAWDS 16 exhibits ab- 
sorption features of He I, but shows no evidence for 
hydrogen. We therefore classify this object as a DB 
WD. As of yet, we have not developed code for fitting 
model atmospheres to DB spectra. We estimate the WD 
temperature based on the equivalent width (EW) of He 
4471 A, as calculated by iKoesteij (|1980D . The EW of 
the He 447lA line (measured from 4350A to 4600A) is 
« ISA; this translates to a temperature of ~ 16000 K for 
log5 = 8. 

It is also possible to use photometry to determine if 
this object is a potential cluster member. Assuming that 
NGC 6633:LAWDS 16 is at the cluster distance and red- 
dening, its photometry (see Table [3|) implies My = 11.56 
and {B — V)q = —0.025. In Table[8l atmospheric parame- 
ters for DB WDs of various masses with My = 11.56 are 
calculated from data presented in IHolberg fc BergeronI 
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TABLE 7 

Cluster White Dwarf Initial Masses 



White Dwarf 








Random Error 


Systematic Error'' 




(Me) 


(Mo) 


(Mq) 





NGC 6633:LAWDS 4 


0.79 


0.01 


2.93 


0.01 


+0.16 
-0.14 


NGC 6633:LAWDS 7 


0.87 


0.01 


3.20 


0.02 


+0.23 
-0.20 


NGC 6633:LAWDS 27 


0.77 


0.01 


3.30 


0.03 


+0.27 
-0.22 


NGC 7063:LAWDS 1 


0.37 


0.01 


8.17 


0.16 


+ <X) 

-1.99 



Prom Table [Gf From uncertainty in cluster age 



((200l . Fr om the table, it is evident that the observed 
photometry of NGC 6633:LAWDS 16 is consistent with a 
cluster DB WD of mass ~ O.6M0 - 0.8Mq, with T^s con- 
sistent with that derived from the EW measurement. In 
addition, Figure [6] shows the DB lying just above the 
O.8M0 DB cooling curve. 

We therefore surmise that NGC 6633:LAWDS 16 could 
be a cluster member. The likely progenitor mass of ~ 
3.5M0 (see Table [U is similar to that of the DA cluster 
member, and if the WD mass is toward the upper end 
of the O.6M0-O.8M0 range, it would fall close to the 
empirical initial-final mass relation. 

If a cluster member, this would represent the third 
He-atmosphere WD to be detected in an open star clus- 
ter, along with LP 475-252 and NGC 2168:LAWDS 28 
(jWilliams ct al. 2006). As such, this WD may be an im- 
portant point in understanding the purported difference 



in the ratio of DA to non-DA WDs in open clusters and 
in the field. Further observations of this object, both 
higher S/N spectroscopy and proper motion determina- 
tion, are needed to confirm whether this DB is a cluster 
member. 

We note that another detected WD in the field, 
NGC 6633:LAWDS 14, is also a cool DB WD. This 
object was shown to be foreground to the cluster in 
iReimers fc Koesteil (|1994f ). 



5.2.3. NGC 7063:LAWDS 1 

distance modulus 
M)v = 9.61 



The calculated distance modulus to NGC 
7063:LAWDS 1 is (m - M)v = 9.61 ± 0.04, less 
than 2f7 away from the cluster distance modulus of 
(m — M)v = 9.41, and is therefore identified as a cluster 
member. However, its mass (O.37M0) is very low for 



its initial mass 



)+oo 



-2.0^^0 



We therefore conclude 
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Fig. 8. — Physical parameters and photometric distance moduUi for WDs in the field of NGC 6633. Points with Icr error bars (often 
smaller than the points) indicate our spectral fits. Horizontal, dashed lines are cooling tracks for WDs with masses of O.4M0 (bottom) to 
1.2M0(top) in O.lM0intervals. Nearly- vertical, dotted lines are lines of constant WD cooling ages from 1 Myr (right) to 1 Gyr (left) in 
logarithmic intervals of 0.33 dex. The heavy, nearly-vertical curve indicates where WD cooling ages equal the cluster age: WDs to the left 
of this line are too old to be cluster members. The shaded region indicates the region of T^ff, logg space where cluster white dwarfs would 
have photometric indices consistent at the 2a level with the observed WD's photometry. This region includes errors in the cluster distance 
modulus, foreground extinction and the observed V-band photometric error. If a WD lies in the shaded region, it is photometrically 
consistent with cluster membership; if it lies outside this region, it is inconsistent with cluster membership. 
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that, if a cluster member, NGC 7063:LAWDS 1 is likely 
a He -core WD and the product of binary evolution 
(e.g.. Ilben &: Liviolll993D . If this hypothesis is correct, 
then follow-up o bservations inay be able to detect the 
companion (e.g.. iMarsh et aLlllOOSl ) or evidence for an 
unseen companion, such as radial velocity variations in 
the core of Ha. This system could then be useful for 
constraining various time scales in common-envelope 
evolution. 

5.2.4. NGC 7063:LAWDS 6 

Our spectral fits to NGC 7063: LAWDS 6 give T^s = 
11,100K and logt; = 7.00, results typical of attempt- 
ing to fit A stars with our routine. However, qualitative 
analysis of the spectrum suggests that the surface grav- 
ity is higher than log (7 = 7. This discrepancy is due to 
a nearby bright star in the slit hampering the sky sub- 
traction to such a degree that a reliable fit cannot be 
obtained. 

If NGC 7063:LAWDS 6 is a cluster member, it would 
have Mv = 12.5 and {B-V)q = 0.33. The range of effec- 
tive temperatures for a WD of this My goes from 8000K 
(for MwD = OAMq) to 13,000K (for Mwd = I-OMq), 
with corresponding cooling ages of 700 Myr and 834 Myr, 
respectively. If the WD is very high mass (1.2Mo), then 
this Mv corresponds to T^s = 20,000K and TcooI = 
511Myr. The color favors lower Teff(~ 7250K). It there- 
fore appears likely that this WD is older than the star 
cluster and so not a cluster member. 

5.3. Expected number of cluster white dwarfs 

The number of expected cluster WDs can be estimated 
by normalizing an assumed IMF to the number of known 
cluster stars. Bas ed on the proper motion study of NGC 
6633 by Sanderi (jl973[ ). the cluster has 45 likely mem- 
bers (P > 50%) brighter than V = 12, corresponding to 
My ~ 3.5. By assuming a Salpeter IMF, an upper mass 
limit to WD progenitors of BM^ n , and by using th e Monte 
Carlo calculations described in iWilliamsl (|2004f ). we cal- 
culate that 7.8 ± 3.0 WDs should be detectable in NGC 
6633, one of which should be binary. This compares with 
< 4 cluster WDs found in this paper. The same calcu- 
lation for NGC 7063, which has 17 likely member stars 
brighter than My w 4.7 in the portion of the field im- 
aged by PFCam, results in 1.1 ± 1.1 cluster WD being 
detectable in the field. 

Thus, like the Hyades, NGC 6633 shows a likely deficit 
of WDs that cannot be explained by WDs being hidden 
in binaries. The NGC 6633 imaging contains almost the 
entire cluster; the CFH12 K field c o vers 4 2' x 28' centered 
on the cluster core, and iSanderd (|1973D detects proper 
motion cluster members over a diameter of ~ 40'. How- 
ever, it is possible that some number of WDs lie out- 
side our imaged area. This areal coverage issue has been 
found to be the li kely cause of the appare nt deficit of 
WDs in Praesepe (|Dobbie et al.ll2004l . l2006h . 

5.4. The initial-final mass relation 

If we assume that both candidate cluster DA WDs are 
indeed members of their respective clusters, we can place 
these two objects on the empirical IFMR. In Figure [TOj 
we plot these two point s along with the data collected 
bv lFerrario eFall (|2005l) . NGC 6633:LAWDS 27 is seen 



to faU on the iFerrario et"al] ()2005D relation, while NGC 
7063:LAWDS 1 falls well off the relation. This can be 
readily explained if NG C 7063 :LAWDS 1 is the product 
of binary evolution (see ij5.2.4p or not a true cluster mem- 
ber. NGC 6633 has a lower metallicity than the Hyades 
and Praesepe, yet the NGC 6633 WD occupies the same 
region of the initial-final mass relation as the WDs in the 
Hyades and Praesepe. This suggests that metallicity, at 
least over this small metallicity range, has little impact 
on the WD mass, though the observational scatter is still 
quite large and can easily swamp a small signal. 

5.5. Three DAZ White Dwarfs 

At least three WDs in this sample show potential Ca II 
K absorption: NGC 6633:LAWDS 8 (EW= 0.6 A), NGC 
6633:LAWDS 15 (EW= 0.4A), and NGC 6633:LAWDS 
27 (EW=0.6A). While these features could be interstellar 
in nature, the other WDs in the cluster field at similar 
or gr eater distances have no discernible Ca II absorp- 
tion. IZuckerman et all (|2003f l determine that ~ 25% of 
DA WDs show Ca II K absorption, if the spectra are of 
high enough S/N and sufficiently high resolution. Our 
detection of Ca II in 3 of 11 DA WDs is fully consis- 
tent with that number. We therefore classify these three 
objects as spectral type DAZ. 

The origin of the metals in DAZ WDs remains contro- 
versial, with explan ations including accretion from the 
interstellar mediuni (|Wesemaellll979t) , cometary impacts 
(jAlcqck et al.lll986j), and accretion of asteroidal mate- 
rial (jGraham et al.l[T990l ). However, recent data from 
the Spitzer Space Telescope and the IRTF reveal that 
several D AZs with high [ Ca/H] have circu mstellar debris 
disks (Rea ch et alll2005HKihc et al.ll2006[ ). Two of these 
DAZs, NGC 6633:LAWDS 8 and NGC 6633:LAWDS 27, 
have Toff in the range where a debris disk may be de- 
tected. 



6. CONCLUSIONS 

We have obtained photometric and spectroscopic data 
for white dwarfs (WDs) in two sparse open clusters, NGC 
6633 and NGC 7063. For the open cluster NGC 6633, 
we determine a main-sequence turnoff age of 560j:^^ Myr. 
We select 32 candidate WDs based onB — V photometry. 
Spectroscopic follow-up of 22 of these candidates finds 9 
of these to be WDs, 7 with spectral types DA or DAZ and 
2 of spectral type DB. One DA and one DB have distance 
moduli consistent with cluster membership. Two other 
DAs are « 0.75 mag over-luminous for cluster member 
WDs; these are potential double-degenerate cluster WDs 
which, if truly binary, would each have combined masses 
above the Chandrasekhar mass. The NGC 6633 DB WD 
is the third known cluster WD with a He-dominated at- 
mosphere. 

For the open cluster NGC 7063, we estimate a main- 
sequence age of 125t^^ Myr, assuming solar metallicity. 
Using t/SUphotometry, we identify nine candidate WDs. 
Spectroscopy of 7 candidates confirms 5 WDs, though 
only one is consistent with cluster membership. This 
WD has a mass of O.37M0, and so is likely a He-core 
WD resulting from binary evolution of two massive stars, 
with the WD's progenitor mass M = SAlt^^MQ. 

The NGC 6633 DA WD and both potential double 
degenerates are consistent with the existing empirical 
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Fig. 9. — Same as Figure[8] except for WDs in the cluster NGC 7063. 



TABLE 8 

Photometrically-derived Quantities eor NGC 6633:LAWDS 16 



1.2 - 



Assumed AfwD {B — V)o 




log g Tcool ^'U 






(Mo) 


(K) 


(Myr) (A/q) 


(Mo) 


1 



0.4 


0.059 


11000 


7.660 


322 


3.72 


+ 0.47 
-0.34 




0.5 


0.008 


12215 


7.868 


305 


3.62 


+ 0.41 
-0.31 




0.6 


-0.026 


13341 


8.017 


310 


3.64 


+ 0.43 
-0.32 


o 


0.7 


-0.056 


14656 


8.163 


286 


3.52 


+ 0.36 
-0.28 




0.8 


-0.076 


16293 


8.320 


271 


3.45 


+ 0.33 
-0.26 


S 


0.9 


-0.088 


18469 


8.480 


250 


3.36 


+ 0.29 
-0.24 




1.0 


-0.110 


21973 


8.650 


205 


3.20 


+0.23 
-0.20 




1.1 


-0.158 


30801 


8.830 


102 


2.92 


+0.15 
-0.14 





initial-final mass relation. Most of the points in the 
initial-final mass relation in the region occupied by the 
NGC 6633 WD(s) are from the super-solar metallicity 
open clusters Praesepe and the Hyades, so the agree- 
ment of the sub-solar metallicity NGC 6633 WD(s) with 
these points suggests little metallicity dependence in the 
initial-final mass relation over this metallicity range. 

Three of our 1 1 DA WDs show Ca II absorption lines 
and are therefore classified as spectral type DAZ. The 
ratio of DAZs to DAs is consistent with other published 
studies indicating that ^ 25% of hydrogen-atmosphere 
WDs are DAZs. 

One open issue is the membership of the reported clus- 
ter WDs. While the observed distance modulus is an 
important diagnostic, we cannot rule out that the clus- 
ter WDs are interlopers from the field. An indepen- 
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Fig. 10. — T he empirical in itial-final mass relation. Filled tri- 
angles are from [F^rrario et al.l 112005): the gray error bars do not 
include the uncertainties in cluster ages. The dashed line is the 
linear fit to the IFMR, from that work. The filled squares are the 
two cluster DAs from this work; the open circles are the potential 
double-degenerate binaries. Errors in final mass are smaller than 
the point size and are therefore not shown. The indicated errors 
in initial mass are Iff errors in the cluster ages. 



dent method of determining cluster membership, such 
as proper motion measurements for these WDs, would 
greatly assist in this regard. 
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